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Abstract 

The anharmonic Bloch oscillations of a light beam in the array of optical waveguides are considered. The 
coupling modes model (CMM) with the second order interaction is used to describe the effect analytically. 
The formula obtained predicts explicitly the path of the optical beam, in particular, the positions of the 
turning points are found. A total agreement of this formula with the numerical simulation is confirmed. 
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I. INTRODUCTION 


Today the arrays of coupling optical waveguides attract the increasing interest since they allow 
to control the behavior of the optical signal effectively. Besides, such systems exhibit many inter¬ 


esting phenomena similar to those that may occur 


Zener tunneling 


in solids jl|, 


such as Bloch oscillations 


1CH12|. Anderson localization 


jl3, 14|, dynamic localization jl5, 16^ etc. 


ffl, 


Usually, for the numerical simulations of the optical phenomena in the arrays of interacting 
waveguides the coupling modes model (CMM) is used. This model describes the behavior of light 
quite exactly if the coupling is enough weak. 

But the theoretical research can be considerably simplified if the numerical simulation can 
be replaced with an analytical formula. In particular, for the periodical arrays of the identical 
waveguides the equations of CMM can be solved analytically (l7|. The analytical solutions were 
obtained for one-dimensional and two-dimensional arrays, and the long-range interaction was taken 
into account. 

Another problem which can be solved analytically is the array of waveguides with the optical 
properties (such as radii or refractive indices) gradually changing as one passes from one waveguide 
to another. This problem was considered in 2|. One of the main results of this work is the analytical 
expression for the path of the optical beam. The obtained formula predicts that the path of the 
optical beam takes the periodical oscillating form. This analytical result was confirmed with the 
experiments 3-6]. This phenomenon is known as optical Bloch oscillations. It is the optical analog 
of the Bloch oscillations of an electron in an ordinary crystal placed in an external electric held. 

Generally, the optical Bloch oscillations, as well as the other optical effects, are investigated for 
the plane arrays of the waveguides, since such arrays are quite simple for fabrication. Besides, such 
systems are the most convenient for the numerical simulation, since the coupling of the adjacent 
waveguides only should be taken into account. 

However, today the increasing attention is drawn to the more complex structures, namely zigzag 
arrays of waveguides (see Fig. [[]). Such arrays are interesting due to the significant role of the 


second order coupling. In particular, in works 


18 


20] the influence of the second order coupling 


on the effects of diffraction and the formation of solitons was investigated. Besides, recently the 
investigations of the optical Bloch oscillations in zigzag arrays were published 21, 22|. It was shown 
that in such systems the signal propagates along a complex periodic trajectory. This phenomenon 
is known as a non-trivial, or anharmonic Bloch oscillations. 
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Figure 1: The scheme of the zigzag array of waveguides. 


The anharmonic Bloch oscillations are the subject of our work. In this paper we generalize the 
results of 2] to the case of the zigzag arrays. We consider the system of equations of the CMM 
taking into account the second order coupling. We find the analytical solution of this system of 
equations and obtain the formula for the path of the optical beam. The obtained formula allows 
to predict the geometric parameters of the anharmonic Bloch oscillations, such as the period of 
the optical beam path and the positions of the turning points. We demonstrate that the analytical 
formula is consistent with the numerical calculations. 


II. ANALYTICAL DERIVATION OF THE ANHARMONIC BLOCH OSCILLATIONS 
BY MEANS OF THE COUPLING MODES MODEL 


Here we represent the analytical derivation of the optical beam path in a zigzag array. Our 
derivation is based on the coupling mode equations with second order coupling. Note that the 
calculation below is a generalization of that represented in [hi], where the ordinary Bloch oscillations 
are considered. 

We start from the equation of the coupling modes 

(i-^Z + Ad 0) + aj'j cij (z) + 71 + a j+1 (z)^j + 7 2 (a i _ 2 (A) + a j+2 {z)^ = 0. (1) 


Here 


Pj 0) = AS 0) + a 3, 


( 2 ) 


For the simplicity, we assume that /5 q = 0. Note that / 3 q can be easily excluded from the 


equations by the replacement aj(z) = a'Az ) e J/3 ° >2 . 


— + cij(z) + 71 ^aj-i(z) + aj + i(z)^j + 72+ aj + 2(z)^j — 0. 


( 3 ) 
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We perform a Fourier transform of cij(z): 


a(k, z ) = 






(4) 


aj(z) = 


V2^ 


dk a(k , 0 ) e 


,*fcj 


(5) 


From Eq. (J3]) one can obtain an equation for a(k, z): 

^i -—f ^ a ~^j S(A;, z ) + 2^71 cos k + 72 cos 2k^j a(k, z) 
Let us consider the eigenmode with a certain /3: 

ap(k, z) = dp(k) e^ z . 

Substituting (J7J) to (JSJ), one obtains: 

i^ a ^k ~ ^ COS ^ ^ 2 COS 5g(A;) = 0 . 

Equation (J8j) possesses the following solution: 

ap{k) = exp < — (jik — sin k — 72 sin 2 A: 


= 0. 


( 6 ) 


(7) 


( 8 ) 


(9) 


The function dp{k) has to be periodic in k. Thus, /3 = an, where n is integer. The function 
dp (k) for a certain n is 


a n (k) = exp <( — ink H-— sin k + — sin 2k 

a a 


( 10 ) 


Now let us turn to the general solution of Eq. ©. Any solution a{k,z) can be represented as 
a superposition of eigenmodes (J7J) : 


l {k,z) = Y J Cnd n (k)e ic 


The coefficients C n can be obtained from the boundary conditions a{k, z 
CD it follows that 

d\k) = Y J C n~a n (k) = 


= 0 ) = 


C n exp <1 — ink H-— sin k H—— sin 2k 

a a 


*72 
a 


( 11 ) 

a°(k). From 

( 12 ) 
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Thus, 


C n — [ — a°(k) exp {ink -— sin k — — sin 2k 

2tt la a 


(13) 


Substituting (1T5T) to (fT2j) . one obtains: 

7 r 

a(M = £ f^At:' mik ’ 

n J 

— IT 

| ^ s i n _ s j n _)_ Yll ^ s i n 2k — sin 2k' 


'-k-\-Oiz) 


(14) 


x exp 


Since 


= 2ir S(k’ -k + az), 


(15) 


Eq. (TT4li takes the form: 

f 2*7 i ( . 


x exp 


a(k, z) = a°(k — az) x 
| —— ^sin k — sin (k — az)J + — ^sin 2k — sin 2(k — az)^j 


(16) 


Let us assume that the function a°(k) possesses a sharp peak near some ko. Then, the function 
at the right-hand side of Eq. (ITU can be expanded into a Taylor series around (k — az) — fc 0 : 


a(k, z) = a°(k — az) e 


i<j)(z)+i(k—az—k o) 


(17) 


where 


4>{z) = —— ^sin (ko + az) — sin k 0 j + — ^sin 2 (k 0 + az) — sin 2k^j , ( 18 ) 

ip(z) = —^cos (ko + az) — cos k 0 j H—— ^cos 2(fc 0 + az) — cos 2k 0 j . (19) 

Substitute mm to (J5j): 


aj(z) = e i< t > ( z )+ i ( k o+ az )i 


dk 

\Z2tt 


d°(k - az) e hfc-^-fco)(i+h(d). 


( 20 ) 


The intensity of the optical beam at the j-th waveguide is defined as Ij(z ) = |a.,-(z)| 2 . Thus, 

2 


Ij(z) = 


iT—ko—az 


— TT—kQ—OiZ 


dk\ 

V2^ 


d°(k 0 + h) e iklU+ ^ z)) 


( 21 ) 


where k\ = k — az — ko. 
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Since a°(k) has a sharp peak at k — ho, the dependence of the integral limits on z does not 
affect the integral value. So, we can write 


W*) = I U + i’(z)) 


( 22 ) 


Thus, the trajectory of the optical beam is defined with the equation j + ip(z ) = const. This 
equation results in 


j(z) = j 0 -— ^cos (ko + az) — cos fcoj-— ^cos 2(ho + az) — cos 2. (23) 

Eq. (1231) represents the trajectory of the optical beam in the zigzag array of waveguides. 


III. QUALITATIVE ANALYSIS OF THE ANALYTICAL FORMULA 

Let us analyze the properties of the trajectory described with Eq. (1251) . For the simplicity, we 
consider the special case jo = 0 and ko = 0. For this case Eq. (1221) takes the form 

j(z) = — — ^1 — cos az'j H—— ^1 — cos 2 az^j . ( 24 ) 

It is obvious that j(z) is the periodical function with the period 2 ix/a. Below we consider the form 
of the trajectory in a single period, 0 < z < 2 tt/ a. We assume that the coupling constants 71 
and 72 are of the same sign (both positive or both negative). This assumption is consistent with 
experiments. 

Let us find the turning points of the trajectory. For this purpose we have to solve the equation 

dj ( 472 \ 

— (z) = 27 ! sin az + 472 sin 2 az = 27 ! sin az 1 H - cos az) =0. ( 25 ) 

dz \ 71 J 


It follows from (1251) . that the turning points are defined with two equations: 

sin az = 0 , 

7i 


cos az = 


472 


(26) 

(27) 


We should consider two different cases, namely 71/472 > 1 and 71/472 < 1. 

Let us begin with the case 71/472 > 1. For this case Eq. (j27l) has no solutions, and all the 
turning points are defined with Eq. (1251) only. In a single period 0 < z < 2n/a Eq. (1251) possesses 
two solutions, Zo = 0 and z u 4 = 7 r/a. Thus, in one period of the trajectory there are two turning 
points. The schematic image of the trajectory is represented in Fig. [2](a). In this figure, the 
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Figure 2 : Schematic image of the trajectory described by Eq. ( 1241 ) for two different cases: 

(a) 71/472 > 1 ; (b) 71/472 < 1 . 


turning points are designated with the letters 0 and A. Substituting the values zo and za to Eq. 
(1241) . one can find the number j of a waveguide where the trajectory changes the direction: jo = 0 
and j A = 471 /a. 


Now let us turn to the other case, 71/472 < 1. The schematic image of the trajectory for this case 
is represented in Fig. [2]^b). For this case, the both equations (j26l) and (j271) possess the solutions. 
Thus, in a period 0 < z < 2n/a the trajectory possesses four turning points. Two of them are the 
solutions of (1261) . they are designated as 0 and A, and their coordinates { zq , jo} and {za, Ja} are 


obtained above. Two other turning points are denoted with B and C. Their coordinates z are the 


1 


a 


7i 


472 


solutions of Eq. (J27J): Zb = — ( vr — arccos - and Zc = — ( vr + arccos - . Substituting zb 


a 


7 i 


472 


and zc to (1241) . one obtains js = jc — 


472 ( 71 \ 

In = = - I 1 -|-I 


a 


472 J 


IV. COMPARISON OF THE ANALYTICAL FORMULA WITH THE NUMERICAL 
CALCULATION 


Below we compare the trajectory described with the analytical formula (1231) with the numerical 
calculation based on Eq. ([[]). For this purpose we take the same parameters as in the paper 211. 
The calculations are produced for a = 0.2, 71 = 1, and for two values of the parameter 72 = 0 and 
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Figure 3 : The trajectory of the optical beam. Left: 72 = 0 ; right: 72 = 0 . 7 . 


0.7. 

For the numerical solution of Eq. (JT]) we have to specify some boundary conditions aj(z) for 
z — 0. We take the boundary conditions in the form of the Gaussian beam 


%( 0 ) 


= e 


ti-i n ) 2 


(28) 


where jo = 10 and a = 4. The Fourier transform a°(k) = a(k, z = 0) defined by Eq. (J4| takes the 
form 

a 0 (k)=e~ s ^. (29) 

This function possesses a sharp peak near k — 0. Thus, in the analytical formula (1231) we take 

k 0 = 0. 

The results of our calculations are represented in Fig. [3] The white trail illustrates the Gaussian 
beam path calculated numerically by means of Eq. (JT]). The dotted line is the trajectory calculated 
with the analytical formula (l23lh The comparison of the numerical and analytical calculations 
confirms that the analytical formula describes the trajectory of the Gaussian beam exactly. 
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V. CONCLUSION 


In this paper we investigated the anharmonic Bloch oscillations by means of the coupling modes 
model. For this problem, the equations of the CMM can be solved analytically. This circumstance 
allowed us to obtain the analytical expression for the path of the optical beam. 

However, to use the obtained formula, one has to define the parameters of the CMM, namely, 
propagation constants and coupling constants. Typically, the values of these parameters are ob¬ 
tained experimentally. 

Recently we proposed a method to calculate the parameters of the CMM by means of the 
multiple scattering formalism |. This allows to simplify significantly the investigation of the 
optical Bloch oscillations. 

Besides, let us notice that the CMM is approximate model, and can be used only for the weak 
interaction between the waveguides. Otherwise, one should use the exact methods of numerical 
simulation. In this case, the analytical formula obtained in this work may be inapplicable. 
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